All-optical steering of light via spatial Bloch oscillations in a gas of three-level atoms 
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A standing-wave control field applied to a three-level atomic medium in a planar hollow-core 
photonic crystal waveguide creates periodic variations of linear and nonlinear refractive indexes 
of the medium. This property can be used for efficient steering of light. In this work we study, 
both analytically and numerically, the dynamics of probe optical beams in such structures. By 
properly designing the spatial dependence of the nonlinearity it is possible to induce long-living 
Bloch oscillations of spatial gap solitons, thus providing desirable change in direction of the beam 
propagation without inducing appreciable diffraction. Due to the significant enhancement of the 
nonlinearity, such self-focusing of the probe beam can be reached at extremely weak light intensities. 

PACS numbers: 42.65.Tg, 05.45.Yv, 42.50. Gy 



I. INTRODUCTION 

Beam steering is one of the most important technolo- 
gies in the modern optics due to its numerous applica- 
tions in such fields as optical imaging, laser machining, 
and free space communication. Various physical mech- 
anisms have been explored for deflection of beams by 
inducing refractive index gradient [H. They include me- 
chanical motion Q, thermal gradient the acousto- 
optical interaction Q, and the electro-optic effect @. 
Fast beam steering in photonic crystals [6[ and phased 
arrays Q were also proposed. 

The direction of the light propagation can also be 
changed by another beam of light through interaction 
with matter. In particular, a medium exhibiting elec- 
tromagnetically induced transparency (EIT) Q can pro- 
vide large probe-beam deflection since the refractive in- 
dex changes significantly near the transparency center 
Q . Among the related studies we mention the refractive 
index measurements by probe-beam (Toj . investigation 
of slow light deflection by magneto-optically controlled 
atomic media II ll . and exploring a scheme of all-optical 
beam steering [TJI]. Yet another example is the electro- 
magnetically induced waveguiding which uses the control 
field as a waveguide to confine the probe field (f3j . 

However, the schemes proposed in the most of the pre- 
vious studies are restricted to the linear regime. They 
usually result in a spread of the probe pulse because the 
refractive index gradient depends on both the probe fre- 
quency and the spatial coordinates. Fortunately, such 
spread can be suppressed by the enhancement of the wave 
localization through the formation of solitons. Large in- 
trinsic nonlinearity in an EIT-based media allows for ex- 
istence of probe beam solitons with extremely weak light 
intensities [f4| • More specifically, ultraslow solitons flol . 
spatial solitons 



I7[ at low light 



[16j . and gap solitons 
intensity can exist in such kind of media. 

In this work, we propose a scheme to achieve 



effi- 



cient all-optical steering of light in a resonant three-level 



atomic system under EIT regime. The scheme is based 
on the phenomenon of nonlinear long-living Bloch oscilla- 
tions (BOs) of gap solitons, recently reported in Ref. fl8j |. 
The system at hand is a gas of A-atoms loaded in a pla- 
nar waveguide created by two photonic crystals. The 
control field used in the system is a standing-wave laser 
beam that originates a linear force, as well as linear and 
nonlinear optical lattices (OLs) affecting a weak probe 
beam. OLs in such a system are flexible: their param- 
eters can be adjusted either by changing the geometry 
and/or intensity of the control field, or by varying one- 
and/or two-photon detunings. 

Briefly, the mechanism of the steering consists of two 
ingredients: the spatially periodic linear force induces 
dynamics of the probe beam in the transverse direction, 
while the combined effect of the lattices controls the di- 
rection of the beam propagation introducing desirable 
and controllable deviations. More specifically, the linear 
OL provides the band structure necessary for existence 
of the Bloch states, while the nonlinear lattice is used to 
control the stability properties of the Bloch states neces- 
sary for existence of gap solitons at both edges of each 
band [llj]. Due to the enhancement of the Kerr nonlin- 
earity induced by the EIT mechanism, stable spatial gap 
solitons, representing the beam channels, can be formed 
even subject to extremely weak probe light intensity and 
subsequently strongly deflected without undergoing ap- 
preciable deformations and attenuations. 

The paper is organized as follows. In the next section, 
the model is introduced. In Sec. Ill, the nonlinear equa- 
tion governing the evolution of the probe field amplitude 
is derived. We show that stable spatial gap solitons can 
exist under properly chosen parameters of the control 
field. We also show how one can implement an efficient 
all-optical steering of the gap solitons. All results in this 
work are obtained under a set of experimentally feasible 
parameters. The outcomes are summarized in the Con- 
clusion. 
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FIG. 1: (a) The geometry considered in the present work: 
the direction of propagation of the probe beam (^-direction) 
is orthogonal to the OL axis (a:-direction). 9 C = arctan(£;/£; c ) 
is the half angle between the two input control beams and 
9 P is the output angle of the probe beam, (b) Left panel: 
A-type atomic system interacting with two laser beams. f2 p 
and Q, c are respectively the Rabi frequencies of the probe and 
control fields. Ai and A2 are respectively the two-photon and 
the one-photon detunings. Right panel: The waveguide cross 
section. The dotted region denotes the core filled in with the 
atomic gas and limited by the Bragg mirrors. 26, h a , and ht 
are the thicknesses of the core and different materials of the 
Bragg mirrors, respectively. 



II. THE MODEL 
A. Preliminary arguments 

We consider a cold gas of lifetime-broadened A- 
type atoms, loaded into an anti-relaxation-coated planar 
hollow-core photonic crystal waveguide, schematically 
shown in Fig. [TJ (It is relevant to note that both room- 
temperature and ultracold atoms have been successfully 
loaded into hollow-core photonic crystal fibers and 
have further been used to study EIT [2(| and all-optical 
switching |2l|). As in any A-system, the atoms can pop- 
ulate three states: the ground state |1), the excited state 
1 2), and the low-energy state |3) [see Fig. []] (b)]. The 



transitions between the states |1) and |3) are forbidden. 

For the particular choice of the atomic gas and geom- 
etry of the system, we take into account the two facts. 
First, in order to contribute to the nonlinear polarization 
the matrix elements of the transition between the lower 
states and the excited state must be nonzero, and hence 
the population of the excited state must also be nonzero 
(see e.g. [TtJ )■ Thus, one has to explore a gas with weak 
atomic losses due to spontaneous emission from the ex- 
cited state. A particular system having such properties 
is the laser-cooled alkaline earth metal atoms, such as 
strontium atoms [22[ where the atomic states |1), |2), 
and 1 3) can be chosen as ^o, and 3 Pi, and 3 Pq, respec- 
tively. One of the advantages of such a system is that 



the atoms possess a long life time even in their excited 
states: the life time of the state |2) is about 21.3 fis cor- 
responding to the decay rate 47 kHz. Therefore, in what 
follows we use this system for the sake of the numerical 
studies (see Figs. [2j EH and©. 

The second fact to be taken into account is the small- 
est possible size of the system. Our aim is to execute 
effective steering of light using OLs. In order to make 
the periodicity appreciable in a finite structure, one has 
to require the system to have the width of at least ten lat- 
tice periods or more. Hence, the desirable OL should be 
created by the beam with the highest possible frequency. 
In our case, there are only two electromagnetic waves 
in the system, which provide the transitions lf)2 and 
3 <-> 2, respectively. Thus we are interested in the geom- 
etry where the strong control field connects the ground 
|1) and the excited |2) states, while the probe weak wave 
links the states |1) and |3). In the case of the stron- 
tium atoms [22| the wavelength of the probe laser beam 
is X p = 52 /xm while the length of the coupling field is 
A c = 689 nm. Requiring the modulation of the nonlinear 
polarization to have at least 20 lattice periods, we con- 
clude that the width of the waveguide should be of order 
of 14 jjxa.. 



B. The geometry of the system 

Now we can describe the geometry of the system at 
hand as shown in Fig. [T] The waveguide created by 
two parallel Bragg mirrors is placed in the (x, z)-plane. 
A strong control field consisting of a superposition of 
counter-propagating waves with the frequency u> c and the 
wavevectors k c ± k = (0, 0, k c ) ± (fc, 0, 0) (k <C k c ), cou- 
ples the excited state |2) and the ground state |1). The 
excited state |2) is also coupled to the low energy state 
1 3) by a weak probe field of the frequency u> p and the 
wavevector k p = (0, 0, k p ) (i. e. it propagates along 
z-axis). We consider a planar waveguide with the trans- 
verse (i.e. along y-axis) width small enough so that only 
one transverse mode of the probe beam needs to be con- 
sidered (see Appendix [5} ■ We concentrate on the 0— th 
mode designated as so(y)- Therefore the probe field can 
be presented in the form 



E p (r,t) = e p £ p {x lZ )s {y)e l{kp2 



-Ul p t) 



c.c 



(1) 



where, e p is the polarization vector and £ p {x,z) is the 
slowly varying envelope of the probe field. 

The frequency (amplitude) of the control field w c 
(E c (r,t)) is much higher than that of the guided probe 
beam u> p (E p (r,t)). This leads to a the effective refrac- 
tive index of the Bragg mirrors for the control field in the 
claddings to be of order one (see Appendix |A"|) . allowing 
us neglecting its transverse distribution (we however em- 
phasize that this design is introduced only for the sake 
of simplicity). The control field can be presented in the 
form 



E c (r,t) = 2e c £ c (x,z)f(x) cos{kx)e^ z -^ 



(2) 
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Here, e c stands for the polarization vector, the function 
f(x) describes the slowly varying distribution of the con- 
trol field in the x direction, and £ c (x,z) is the complex 
amplitude of the control field. In experiments, the de- 
sired space dependence of f(x) can be achieved with the 
help of beam masks. 

The total electric-field can be written down as 
E(r, t) = E p (r,t) + E c (r,t) and considered classically. 
Large difference between the frequencies of the probe and 
control fields (ui c 3> Wp) allows us to split the equations 
for E p and E c , while large difference in the field ampli- 
tudes (E c E p ) allows us to consider E c to be constant 
in the equation for the probe beam . Now the equation 
governing E p can be written down as follows: 



rf 1 d 2 
V 

c 2 Q t 2 



E p = 



1 d 2 



e c- 



dt 2 



P. 



(3) 



Here P is the probe beam polarization defined by the 
properties of the atomic gas in the waveguide core and 
by the properties of the Bragg mirrors in the waveguide 
claddings, i.e. 



Here f2 p = p 2 3,£ p /h and il c = p2i£ c /K are the Rabi fre- 
quencies of the respective fields A2 = (0J2 — ^3) — w p , and 
Ai = (cl>3 — Wi) — (uj c — u)p) present one- and two-photon 
detunings, respectively (see Fig. [IJ. In Eqs. ©, we keep 
only the dissipation of the state |2) due to spontaneous 
emission, by adding the decay rate 72, phenomenologi- 
cally. Usually, the decay rates of the lower states |1) and 
1 3) are much smaller than that of the exited state and can 
be safely neglected. For a particular choice of the atomic 
gas and assuming that the decay rate of the excited state 
is much less than the two-photon detuning (72 -C Ai), 
one may consider the loss of the atoms from the system 
as a small perturbation (see also the discussion in [I3|)- 
Since the two-photon detuning is nonzero (i.e. Ai ^ 
0), the pure dark state (corresponding to %p 2 = 0) can 
not exist. However, we explore the situation where the 
atomic system is close enough to the dark state. To 
this end we look for a solution of Eqs. ([5]) in the form 



i)j(r,t) 



i\t. 



, (r), where the exponent A is slowly vary- 



ing in space and can be obtained as a root of the respec- 
tive characteristic equation. We are interested in a root 
having the smallest imaginary part [24j , i.e. 



P 




(y <-b,y> b) 
(-b<y< b) 



(4) 



Following the standard procedure [23| P c i ac j can be deter- 
mined by calculating the effective refractive index of the 
Bragg mirrors n from Eq. (|A3[) : P c i a d = £o(™ 2 — 1)-Ep. 



C. Nonlinear Polarization 

In order to compute P corc one introduces the bosonic 
field operators tpj, of the states \j) (j = 1,2,3), as 
well as the electric dipole matrix elements pij associated 
with the transitions between the states \i) and \j). Then 
the polarization can be obtained as 



(5) 



The steering of the light we are interested in is per- 
formed through the proper spatial and/or temporal mod- 
ulations of Pcorc- In order to link it directly to the probe 
field E p we have to address the Heiscnbcrg equations for 
the operators ipj describing the atomic medium. This 
can be done neglecting the atomic kinetic energy, what 
leads to the system as follows (its derivation repeats the 
steps outlined say in [l7j |) 

d - 

i-^-tpi = -Ai-01 - 2n*f(x) cos(fca;)-02, (6a) 
at 

d - 

i-^:4>2 = (A 2 - 172)^2 - 2Q c f(x) cos(fca;)-0i 
ot 



.8 r 



-ftpips, 

-Kit- 



(6b) 
(6c) 



(I Ai| < |Ar|) with 



A = A,. + iXi, 



Ai|fip| 



(7) 



AiA 2 + A\n c \ 2 f 2 { x ) cos 2 {k x ) + \n p \ 2 ' 

72 x2 
{ 2, 



A ' — 10 \2 A r 



Then the expectation value (^3^2) can be readily com 
puted from © and ([7]): 



t t x A|^i-A| 2 » 

sV2) (|A| 2 + |fi| 2 )|^ - A| 2 + 4|A| 2 / 2 (x) cos 2 (kx) ' 



Here we introduced the dimensionless functions O 
Q p /Q c and 

A = -^-=A r + ,A i , 



(8) 



(9) 



where 



A,. 



A,- = ^A 2 



6162 +4/ 2 (x) cos 2 (kx) + \n\ 2: 
^A 2 

n\ 2 r: 



as well as the parameters 5j = Aj/|n c | (j = 1,2) and 
72 = 72/I |- Since the probe field (decay rate of the 
excited state) is assumed to be weak in comparison with 
the control field, we have that < 1 (72 < |J1| 2 < 1, 
see Eq. (jTTJ) and the related discussions). 

An important conclusion follows from Eq. (|6cp . the 
quasi-dark state, characterized by the small population 
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of the excited state |2), requires |A| <C |O p | and hence 
|A| « 1. 

In what follows we will be interested in a particular 
choice f(x) = 1 — ax where a is the mask parameter 
considered to be small. More specifically, we will con- 
sider the beam propagation (i.e. distributions of in 
space) characterized by finite deviations from the axis x. 
As we have already announced in the Introduction, the 
mentioned deviation will occur due to BOs of the beam. 
Thus if we assume that the amplitude of BOs measured 
in units fc" 1 [sec Eq. (Tl3|) below] is X, the smallness of a 
is determined by the requirement that aX/k ~ Q 2 <C 1. 
This allows us to achieve further simplification for the 
polarization by defining 



with 



v(kx) = 6162 +4cos 2 (fca;). 



(10) 



Performing the expansion of the right hand side of Eq. 
© in the Taylor series with respect to Q we get: 



Sin 

v(kx) 



8ax cos 2 (kx) 



v(kx) 



8\ - M2 + 2v{kx) 
v{kx) 2 



v(kxy 



n. (11) 



Here we have neglected the H 5 -order terms and assumed 
that in the domain of BOs the condition |S!| 2 <C v(x) is 
satisfied. Notice that the requirement v(x) 7^ holds for 
in the whole space provided 

8i6 2 > or 8182 < -4. (12) 

This condition will be imposed in what follows. 



U 



Aflp 3 2l 2 fe^ 

Henk 2 Q r 



and a' — a/k. The constant Co depending on the par- 
ticular choice of the Bragg mirrors is computed in the 
Appendix [A] 

The obtained dimensionless Eq. (fT3]) implies that all 
the terms are of the order one or less. More specifically, 
it is necessary to require that Uq ~ 1. This last con- 
strain can be satisfied by different means. Notice that 
the parameter k p /k can be experimentally controlled by 
the geometry of the laser beams (see Fig.[T]), we thus fix 
k p /k = 0.2 and adjust other parameters (in particular S± 
and O c ) to satisfy the constrain. For a typical set of data 
explored below, i.e. the atomic concentration Af « 10 14 
cm" 3 and Ai = O c = 1.0 x 10 7 s~ x (81 = 1), we ob- 
tain Uq = 4.7 satisfying the desired order of magnitude. 
Meanwhile, the dissipation is very small (r(£) <C 1) in 
the considered system since 72 = 4.7 - 10~ 3 <C 1. 



III. SPATIAL GAP SOLITONS AND 
ALL-OPTICAL STEERING 



A. Spatial gap solitons 

We study the situation when the effects of the lin- 
ear force, nonlincarity, and dissipation are small enough. 
Therefore we first notice that for F = 0, G = 0, and 
r = one arrives at the familiar framework of the band 
theory, where by means of substitution O oc e -^ K a( 
Eq. (fl3|) is reduced to the eigenvalue problem 



D. The equation for the beam envelope 

We are particularly interested in the beam dynamics 
resulting from the interplay of diffraction and nonlinear- 
ity when the probe field passes through the atomic cloud. 
Therefore, we concentrate on stationary solutions, requir- 
ing fl to be independent on time. Substituting Eqs. (fTJ) 
and (jlip into Eq. ([3]), taking into account (|A1|) . and 
leaving only the leading order terms, we arrive at the di- 
mensionless equation for the slowly varying function O: 



.90 d 2 fl TT ,„ n 
+iT(t)n = 0. 



■tF(on-G(0\n\ 2 n 



(13) 



Here £ = kx and £ = (k 2 /2k p )z are the dimensionless 
independent variables, v(^) = v{kx) is given by (|10[) . and 
we defined 

m = Uo m = Sa'U 0COS ^) 



v 2 (0 



G(e) = C U [8 2 - 8jM + 2i;(fl] r <5 l72 C/ 



«(0 S 



£u a ,q(£) + K a {q)u a , q (£) = 0. 



(14) 



Here C = —d 2 /d^ 2 — J7(£) and M Q , g (^) is a Bloch func- 
tion with indexes a and q standing respectively for the 
band index and for the wave vector in the first Brillouin 
zone (BZ), i.e. for q € [—1,1] in the case at hand. The 
physical meaning of the eigenvalue K a (q) is the propa- 
gation constant. In what follows we concentrate on the 
lowest energy band (a = 0) and therefore, for the sake of 
brevity, omit the band index. 

The effect of the periodic linear force originated by the 
control beam can be described by the equations 



tan6> = 



dK a (Q) 
dQ ' 



dQ 
d( 



-F(E) 



(15) 



where 3 and Q denote the coordinates of the center of 
the Bloch wave packet in the real and reciprocal spaces, 
respectively. 9 is the refraction angle with tan# = 
2(k p /k) tan# p (see Figs. Q] (a) and [3] (a)). A peculiarity 
of our case is that F(£) is an oscillating function, which 
has a nonzero mean value 



1 



(F) = - F(0 = 



4a' V a 



|£i 62 + 4| y /(6 1 5 2 + 4)6162 
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and hence results in the oscillation behavior of the Bloch 
wave packet in the real and reciprocal spaces (i.e. in 
BOs). The dimensionlcss amplitude (along the x axis) 
and period (along the z axis) of such oscillations can 
be calculated as X = AK/2(F), with AK denoting the 
width of the lowest allowed band of the periodic potential 
f/(£), and Z = 2/\{F)\, respectively. 

In order to describe the effect of the nonlincarity we 
employ the standard multiple-scale expansion assuming 
smallncss of the averaged force, i.e. (F) <C 1. This last 
requirement is equivalent to the condition a« fc, which 
is consistent with the constrains on the BOs amplitude 
discussed in Sec. II C provided X < 1. 

Next we introduce the scaled variables (£j,£j) = 
fj?/ 2 (f;, Q (j = 1, 2, ■ ■ •), where \i is a small parameter es- 
timated by fi ~ (F) <C 1 • The probe field can be written 
in the form of the expansion £1 = Yl'jLi M ^ (£oj Ci)j 
where in the arguments of the functions we have in- 
dicated only the most rapid variables. T(£) = ^ 2 7?.(£o)- 
Substituting this ansatz into Eq. (fT3|) and collecting the 
terms at each order of /i 1 / 2 , we obtain a set of equa- 
tions which can be solved order by order (see e.g. |25j). 
Omitting the details we turn directly to the equation 
for the slowly varying amplitude A(XjCi)i which is de- 
fined through the relation fiW = A(x, (2)e lK{QKo u Q (£ ), 
where x = £i — tan[0(£ 2 )] Ci an d the dependence Q = 
(5(^2) and 6(^2) arc obtained from Eqs. (fT5|) . In the 
third order of the asymptotic expansion, 0(/^ 3 / 2 ), we ob- 
tain 



r)A f) 2 A 

t +V(Q)—-G(Q)\A\ 2 A 
0(2 ox 2 



iK{Q)A = 0, (16) 



where the coefficients of the diffraction T> and the effec- 
tive nonlincarity Q are given respectively by 



1 <PK{Q) 
TZ(Q) = f n(0\u Q (0\ 2 d^. 



G(0\u Q (®\*d£, 



We notice that both T> and Q can be either positive 
or negative depending on the signs of detunings and on 
the position of the center of the wave packet (i.e. on 
£2)- Therefore, in a general case BOs oscillations will be 
accompanied by significant change of the width of the 
beam, due to the diffraction. This undesired effect, how- 
ever can be dramatically reduced, if the conditions for 
sclf-focusing of the beam will be satisfied (and hence will 
compensate the diffraction) at any coordinate £ 2 . These 
requirements (which are also referred to as conditions for 
the gap soliton existence or, alternatively, as conditions 
of the modulational instability of the Bloch waves) , are 
well known and read [Ifl [26[ 



V(Q)G(Q) < 0. 



(17) 



If condition (fT7|) is satisfied for all (or almost all) wave- 
vectors Q, the sclf-focusing of the beam will occur along 
the whole (or almost whole) trajectory of the beam. 
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FIG. 2: (a) The white domains correspond to the choice of 
the detuning, such that (|17|l is valid for 90% of wavevectors in 
the first BZ; in the black domains the requirement (|12p is not 
satisfied and thus the approximation Eq. (|13l) fails; the gray 
domains show the parameter region where our approximation 
is valid, but the conditions (|17|l are not satisfied for more than 
10% of the wavevectors in the BZ. The system parameters are 
given in Sec. II A. (b) The coefficients T>(Q) (dashed line), 
Q{Q) (dash-dotted line), and T>(Q)Q(Q) (solid line) vs Q with 
the parameters Uo = 4.7, Ai = fi c = 1.0 x 10 7 s _1 , and 
A2 = —1.0 x 10 s s _1 [the point A in panel (a)], (c) The same 
as in panel (b) but with A2 = —0.6 x 10 s s" 1 [the point B in 
panel (a)]. The dynamics of the Bloch wave packets for the 
points A and B is shown below in the panels (a) and (b) of 
Fig- El 



Thus, as the third step, for achieving the best perfor- 
mance of the device one has to design the system (geome- 
try of the beams or photon detunings) in a way to ensure 
the condition ([TT]) for all Q. The procedure of perform- 
ing this task was described in the previous works (l8l . [2(| 
and here we report only the final result. More specif- 
ically, in Fig. [5] (a), we show the domains on the plain 
(6\, 82) where the optimal conditions are satisfied (empty 
regions), not satisfied (gray regions), and the parameter 
requirement (fT2|) is violated (black regions). In Fig.[2](b) 
and (c), we plot the curves T>(Q), G(Q), and V{Q)Q{Q) 
versus Q with the almost optimal design [point A in panel 
(a)] and non optimal design [point B in panel (a)] of the 
OLs. We mention that always there exists a small range 
of the wavectors where spatial gap solitons do not exist 
(i.e. where T>(Q)Q(Q) > 0). This range, however, can be 
made enough small [less than 10% of the whole BZ, like 
this happens in the point A, sec Fig.[2](b)], and therefore 
has no appreciable destructive effect on the beam dynam- 
ics [see Fig. [3] (a)] . On the contrary, the point B in panel 
(a) corresponds to the case in which the BOs of the beam 
are rapidly destroyed leading to strong defocussing of the 
outcome beam [see Fig. E](b)]. 
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(b) 




FIG. 3: (color on line) (a) The long-living BOs of the beam 
in the almost optimally designed OLs [point A in Fig. [2] (a)] 
obtained by direct numerical integration of Eq. (|13p . The 
amplitude and the period of the BOs are X = 14.47 (= 49 
fim) and Z = 578.98 (= 787 fim). (b) BOs of the wave packet 
in the non optimally designed OLs [point B in Fig. [2] (a)]. In 
the both cases the initial condition is a stationary gap soliton 
located at £/X = 0. The intensity profiles of the output 
beam (blue, dashed line) vs the input beam (red, solid line) 
corresponding to the cases (a) and (b) are respectively shown 
in the panels (c) and (d) obtained at ( = Z. 



B. Long-living BOs and all-optical steering 

Now we turn to the final step of numerical descriptions 
of the BOs of a spatial soliton and of the ight steering. In 
Fig- El (a), we show the dynamics of the Bloch wave packet 
launched normally in the atomic media with the almost 
optimal design of the OLs [point A in Fig. [2] (a)]. The re- 
sult was obtained by employing numerical integration of 
Eq. (fT3|) . The existence of the long-living BOs of the gap 
soliton, with the amplitude and period well matching the 
theoretical estimates, can be observed. Due to the very 
small dissipation, there is only a tiny decrease of the soli- 
ton intensity due to the absorption of the medium. More 
specifically, one can see the highly concentrated beam in 
the output with the direction of propagation determined 
by the nonzero angle 9. For the sake of comparison, in 
Fig-El(b), we show the dynamics of the same Bloch wave 
packet in the non optimally designed OL [The parameters 
used in this panel are the same with those used in Fig. [5] 
(c)]. In this case, the beam undergoes strong spreading 
out due to the alternating diffractive and self-dcfocussing 
regeems. The intensity profiles of the output beam vs 
the input beam are shown in the panels (c) and (d) . One 
observes that there is no apparent deformation between 
the output beam and input beam in panel (c) while an 
obvious destruction appears in panel (d). 

The maximum input intensity of the probe beam in 
Fig. [3] (a) can be estimated by the formula I max = 
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FIG. 4: tan (9 vs 9 C at ( = Z for the almost optimally designed 
OLs [corresponding to point A in Fig. [2] (a) and panel (a) of 
Fig. 0. 

(c/2)eo|E pmax | 2 . Using the parameters given in Fig. [2] 
(b), we obtain I max « 4.5 fiW cm~ 2 . Thus, the gener- 
ation of the spatial gap soliton in the system at hand 
requires only very low input light intensity. We remark 
that the intensity of a single 500-nm photon per nanosec- 
ond on an area of 1 /im 2 is I p h = 0.04 W cm -2 . This 
estimation shows that our system makes it possible to 
manage the weak beams characterized by single-photon 
wavepackcts. This is drastically different from the spatial 
optical soliton generation in a conventional waveguide 
where an input laser pulse with very high peak power 
~ 10 2 kW is needed in order to bring out a sufficient 
nonlinear effect (27[. 

The existence of long-living BOs of spatial optical gap 
solitons can be used to implement efficient all-optical 
steering of probe beams. Different from all linear sys- 
tems where the spread and attenuation of the probe pulse 
are unavoidable because the refractive index gradient de- 
pends on both light frequency and spatial coordinates, 
the scheme proposed here can greatly improve the de- 
vice performance due to the self-focusing of the beam. 
In experiments, the optical masks and atomic cells arc 
usually chosen in advance. However, one can exploit the 
possibility of changing the geometry or parameters of the 
control field to steer the output probe beam. As an ex- 
ample, by changing the angle between the input control 
fields 9 C = arctan(fc/fc c ) (see Fig. [1]), one can efficiently 
control the refraction angle of the output probe beam. 
In Fig. 01 we show the tangent of the refraction angle 
tan 9 versus 9 C with a particular length of atomic cells 
(i.e. C = Z). A wide range of angle, 9 p e [-0.98, 1.19], 
can be achieved with the given parameters. 



IV. CONCLUSION 

In conclusion, we studied (both analytically and nu- 
merically) the dynamics of weak probe optical beams in 
a gas of three-level A-atoms subjected to a control laser 
beam. The standing wave originates periodic modula- 
tions of the linear and nonlinear parts of the dielectric 
permittivity of the medium, which is "seen" by the probe 
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beam. If the amplitude of the control field is smoothly 
modulated, the probe beam undergoes the Bloch oscilla- 
tions. These oscillations correspond to periodic change 
of the direction of the propagation of the beam. We have 
shown that it is possible to choose the parameters of the 
atomic system, say Rabi frequency of the control field 
and photon detunings, to reduce the diffraction of the 
probe beam, and thus the output beam is approximately 
the same with the input one. In this case, the probe field 
distributions is nothing but a spatial gap soliton undergo- 
ing long-living Bloch oscillations. Due to the significant 
enhancement of nonlincarity, such solitons can be formed 
with extremely weak light intensity, below single-photon 
wavepacket level. Thus, through changing the geometry 
or parameters of the system one can control the direction 
of the output beam to carry out an efficient all-optical 
steering of light. 



Acknowledgments 

The work of C.H. was supported by the Fundagao 
para a Ciencia e a Tecnologia (FCT) under Grant No. 
SFRH/BPD/36385/2007 and Estunulo a Invcstigacao 
2009 de Fundagao Calouste Gulbcnkian. The research 
of VVK was partially supported by the grant PIIF-GA- 
2009-236099 (NOMATOS). 



Appendix A: The transverse profiles of the probe 
field 



The transverse profiles of the probe field can be deter- 
mined from the eigenvalue problem [23TJ 



^ + ^- k ?) s i(y)=°A\y\<b) (Ala) 
^ + ^ - kf) si(y) = 0, (|y| > b) (Alb) 



tan 

mode index (I = 0,1,...), si(y) is a normalized, i.e 



where uj stands for either u> p or w c , I stands for the 



I^oo dy\si{y)\ 2 = 1, profile of the respective mode, and 
2b is the width of the waveguide core. 

The eigenvalue ki and effective refractive index of the 
Bragg mirrors n can be obtained from the dispersion re- 
lations of the planar waveguide 



tan 



(Z + 1)tt 



2^,2 



kfc 



r, (A2) 



and the expression valid the finite photonic struc- 
tures m 



n(ui) 



c 



(fit — - lnT(w) 



(A3) 



where (fit(ui) and T(w) are the total phase and the trans- 
mittancc of the complex transmission coefficient for the 
structure t(u) = x{oS) + iy{oS) = VTe 1 ^* with <j> t = 
arctan(y/x) ± mir (the integer m is uniquely defined as- 
suming <fi t is a monotonically increasing function) and 
T(u) = x 2 + y 2 . h is the total thickness of the A^-period 
crystal h = N{h a + hb) with h a and hb being the thick- 
nesses of two different materials [see Fig. [T] (b)]. 

As an example we consider the probe field constrained 
of quarter- wave Bragg mirrors, n a h a = n^hf, = A p /4 with 
n a = 1.0 and rib = 1-4 being the refractive indexes of the 
two different materials, for claddings, N — 2, and 2b = 
0.1 mm, we obtain that fco = 1081 cm -1 and n m 0.87. 
At the same time, for the control field we have fco ~ 91144 
enr 1 and the effective refractive index of the cladding is 
n w 1.0. 

The assumption that only s takes into account re- 
quires that the energy for exciting si [corresponding to 
the eigenvalue of Eq. (|Ala.[) ] is higher than that of the 
self-action (nonlinearity) in Eq. (fT3|). With the param- 
eters given in Fig. we obtain G(^)\Q\ 2 Q < (ujp/c 2 — 
k 2 )/k 2 (ki = 873 cm -1 for the probe field), hence we 
can safely neglect the higher transverse modes (i.e. si 
for I > 0). Here, the integral C = \s \ 4 dr] = 73.9. 
However, the much higher frequency of the control field 
allows us to neglect the difference among transverse pro- 
files of the control field in the waveguide core. 
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